Increasing environmental pollution, which causes global warming and endangers human health and environmental degradation, has led to the concern of many supply chain managers and designers. The purpose of this research is to provide a mathematical model for designing the purchase, production, and distribution in a multi-level and multi-product supply chain network such that the environmental impact and total costs of supply chain is minimized and the customers' satisfaction level is maximized. Due to the unspeci, ed demands level, demand uncertainty has been considered in the problem. Regarding the complexity of the proposed mathematical model and di, culties in solving the problem with exact methods in large size, a NSGA II has been proposed. To evaluate the proposed NSGA II, , ve sample instances are generated in di , erent size and solved by Epsilon constraints method and NSGA II. According to the results, the proposed NSGAII is a reliable method to , nd e, cient Pareto frontiers in a reasonable time.